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A hydrostatic weighing system is described that provides independent values of apparent 
mass with a standard deviation of about 0.2 microgram. The definition of ' 'independent" 
used here includes the requirement that the sample under observation be removed from the 
water, dried, and recleaned between measurements. This precision is between one and two 
orders better than existing hydrostatic balances, and permits high-quality density measure- 
ments without recourse to large sample sizes. The system is used as a comparator, hence 
the absolute accuracy of results can be no better than that of the standards used. Data 
are presented which were taken during experiments on two 2-gram pieces of single crystal 
silicon, and the standard deviation in density of a single determination was computed to be 
4.4X10" 7 g/cm 3 (0.22 ppm). 

The system is particularly well suited to detection of slight density changes in small 
size samples. In this service, it is not subject to the requirement for high accuracy standards. 
Data are reported on density changes of about 50 ppm that occurred in a 250 milligram 
ceramic crystal. Two independent determinations on this change differed by only 7 percent. 



1. Introduction 

The density of an irregularly shaped object is 
usually determined by a combination of air and 
water weighing in which the apparent mass of the 
object is determined in both media. The apparent 
mass in air, A, is determined on a balance in the 
usual manner. The apparent mass in water, B, is 
determined by suspending the object from one pan 
of the balance by a fine wire into a body of water. 
The mass, M, and volume, V (the ratio of which is 
density), result immediately from a simultaneous 
solution of 

A = M— p A V (air weighing) 
B=M—p w V (water weighing) 

where p A is the density of the air displaced by the 
object during the air weighing, and p w the density 
of the water displaced by the object during the 
hydrostatic weighing. 

The main limitation in the above scheme occurs 
in the determination of B. There are two reasons 
for this: 

(1) Although we know the density of a particular 
sample of water to almost two orders better than 
we know the density of a particular sample of air, 
water density is about three orders greater than that 
of air. Hence, the uncertainty in p w V is greater 
than that of p A V. 

(2) In water weighing, the variability of non- 
buoyant forces is of much greater magnitude than 
those occurring in air weighing. First, there are 
variable momentum forces associated with turbu- 
lence and convective currents in water which act on 
the object being weighed. Such forces also exist in 
air weighing but to a much smaller extent. Second, 
there is variability in surface tension forces associated 



with the meniscus surrounding the point where the 
support wire penetrates the water surface. Under 
nonideal conditions, these forces may vary by 10 /xg 
or more. 

Thus, the accuracy in a measured value of B is 
limited by our knowledge of p w , and the reproduci- 
bility of the measured value is limited by the 
variability of the nonbuoyant forces mentioned 
above. 

For a given hydrostatic weighing system the ab- 
solute value of the variability (due to nonbuoyant 
forces) is largely independent of load. Hence, to 
achieve improved values of percentage reproducibility 
in B, we must use samples of very large sizes. 

Unfortunately, objects for which density deter- 
minations have the maximum scientific value are 
rarely available in samples of unlimited sizes. It is 
seldom that samples above 2 g are obtainable, and 
they are usually of fractional-gram size. Such 
items as uniformly irradiated crystals, samples of 
maximum purity substances, and balance sensitivity 
weights, usually weigh 250 mg or less. For a 250- 
mg object, a variability of 10 /xg in the nonbuoyant 
forces would result in a variability in B of about 
40:10 6 — a highly unsatisfactory measurement for 
many technical purposes. 

In the Cartesian diver system to be described, it 
is possible to determine B of a 250-mg object to a 
reproducibility (standard deviation) of about 1:10 6 . 
The reason for this is that the two major sources of 
variability in ordinary hydrostatic weighing do not 
exist in the diver system or are minimized: 

(A) No portion of the diver structure penetrates 
a water surface, so uncertainties caused by surface 
tension do not exist. 

(B) The diver system is operated at a nominal 
temperature of 4 °C; hence, nonideal thermosta ting- 
does not result in variability of water density. In 
the absence of density gradients (other than those 
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caused by gravity) convection currents cannot be- 
come established, and turbulence damps out quickly. 
Thus momentum forces are minimized. It is gen- 
erally considered impractical to perform ordinary 
hydrostatic weighing at 4 °C. 

2. Description of the System 

The Cartesian diver is not a new device. A re- 
view of prior art is given by Chiarulli and Chilton 
[l]. 1 Recent work by Haller and Calcamuggio [2] 
and Spaepen [3] has raised the state of the art 
greatly. 

Basically the diver, figure 1, is a hollow, compres- 
sible body with a closed surface to which are attached 
various hooks, guides, and reticles to facilitate load- 
ing, manipulation, and observation. When such a 
body is immersed in water and the pressure in the 
water is increased, the positive buoyancy of the 
body will decrease or increase respectively, depend- 
ing upon whether the compressibility of the diver is 
greater than or less than that of water. If diver 
compressibility is equal to that of water, buoyancy 
is independent of pressure. 



Figures in brackets indicate the literature references at the end of this paper. 



OBSERVED TIP 



UPPER PORTION 
OF BODY IS 
CYLINDRICAL 




LOWER PORTION OF BODY 
PARTIALLY FLATTENED TO 
INCREASE COMPRESSIBILITY 



THIN WALL HOLLOW 
^QUARTZ BODY 



TARE WEIGHT HOOKS 



LOWER FRAME MADE 
OF I mm FIBER 




LOAD BASKET HOOK 



ALINEMENT TIP 



Figure 1. 



Our most useful diver is made of fused quartz 
(vitrified silica) and weighs about 5 g. It displaces 
about 10 cm 3 ; hence in water it has a positive buoy- 
ancy of about 5 g. It has a very thin wall so it is 
quite compressible. In practice it is loaded by the 
object to be measured and enough tare weight to 
cause it to sink slowly in water at 4 °C. The pres- 
sure system, figure 2, is then closed and the water 
pressure reduced slightly. This causes the diver 
volume to expand, thereby increasing diver buoy- 
ancy. This causes the loaded diver to stop falling 
and begin to rise. 

There is a pressure gradient in the water, due to 
the gravity field, so that as the diver rises, it moves 
into a region of lower pressure. This causes addi- 
tional expansion and an increase in buoyance, so the 
diver is accelerated upward. Conversely, when the 
water pressure is increased, the diver body is com- 
pressed which results in smaller buoyance. The 
diver will then sink. As it sinks into a region of 
higher pressure, additional buoyancy is lost and the 
diver is accelerated downward. Thus, the gravity 
pressure gradient in the fluid always accelerates 
diver motion, and as a result, there is no position 
of stable equilibrium. 

On the other hand, had the diver compressibility 
been less than that of water, the reverse set of 
circumstances would have held. Under such con- 
ditions, the pressure gradient in the water would 
always impede the diver motion and a position of 
stable equilibrium would exist. It is possible to 
operate a diver system under either condition, but 
practical considerations such as load range and 
sensitivity led to our choice of a diver more com- 
pressible than water. 

Diver motion is observed through a rigidly 
mounted microscope and a system of windows 
and mirrors. The optical center of this system is 
carefully adjusted horizontally, and the horizontal 
plane within the diver chamber containing the 
optical center is called the reference plane. Although 
(as explained above) there is no level of stable equi- 
librium for the diver in the system, there is a pressure 
which will result in the loaded diver having a density 
precisely equal to the density of the water it displaces 
when the observed upper tip of the diver is in the 
reference plane. This pressure is called the equi- 
librium pressure, P E , and it can be determined by 
averaging the two pressures which cause the diver 
to move upward and downward across a vertical 
range defined by an eyepiece reticle and extending 
about 0.013 cm above and below the reference plane 
in equal intervals of time (usually about 30 sec). 

When the pressure is very close to P E) we are 
justified in assuming that the average velocity, y, of 
the diver, during its transit of the range, is nearly 
linear with pressure. Thus, when the average 
rising velocity, y R , does not precisely equal the 
average falling velocity, y F , an acceptable value of 
P E can be calculated from the formula 

P e P r Vr 



Pf p r vf yR 
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Figure 2. 



Alternatively, there is a quick graphic solution, 
which we used, illustrated in figure 3. 



3. Pressure System 

Our pressure regulator system is a modification of 
the two-bellows barostat developed by D.P. Johnson 
[4] for use in the NBS barometer calibration service. 
It consists of two stainless steel bellows mounted 
one above the other. The upper bellows is evacu- 
ated, and the lower is water-filled and connected 
to the diver chamber. The upper end of the upper 
bellows and the lower end of the lower bellows are 
rigidly attached to the frame. The two bellows are 
connected mechanically together by a freely floating 
separator plate. Under such circumstances, the 
pressure developed in the water in the lower chamber 
is directly proportional to the deadweight loading 
of the separator plate. The effective areas of the 
bellows are such that a 76 kg deadweight load on the 
separator plate results in about 1-atm pressure 
being developed in the lower chamber. The loading 
of the separator plate is conveniently changed in 
increments of about 100 g, from a minimum of 2 
kg (the tare load) to 85 kg, by a semiautomatic 
loading machine developed by David Sklar, Design 
Development Laboratories, for the purpose. A 
small tray on top of the assembly is used to receive 
additional small hand-loaded weights of between 1 
and 100 g (see fig. 2). 

The separator plate is electrically connected as 
the moving plate of a two-plate variable capacitor. 
The fixed plate is rigidly attached inside of the 
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evacuated bellows, about 0.013 cm above the mid- 
point of the 0.025 cm free travel allowed the separator 
plate by its stops. The capacitance between the 
fixed and moving plates is fed into a Wien bridge 
oscillator and discriminator, which provides servo 
power to drive a motor-driven bellows connected in 
the hydraulic system, in such a direction that the 
separator plate is maintained at the midpoint of 
its allowed travel, regardless of loading. 

With this arrangement the two-bellows assembly 
is used as an error signal generator and actual pres- 
sure changes are made by the servobellows. Inas- 
much as the separator plate remains at the same 
vertical level, regardless of pressure, the effects of 
the elastic properties of the stainless steel bellows 
are largely eliminated. 

The elastic properties of one of the divers are 
such that the diver buoyancy changes by 1 p.g when 
the separator plate load is changed by 12 g. Inas- 
much as diver buoyancy on differences is constant 
to about 2/10 pg, it can be inferred that the short 
term stability of drift of the pressure regulator is 
about 2X10 -5 atm/hr. This is consistent with 
other estimates of stability of similar devices [5]. 

As use of the diver system does not require absolute 
knowledge of pressure, no effort has been made to 
calibrate the system in absolute terms. In use, 
pressures developed are stated in terms of the 
kilogram load on the separator plate, which is 
acceptably linear with water pressure. 



4. Use of the Diver — Density Measurements 

If 6?i is the mass of water displaced by the Carte- 
sian diver body at some pressure, then 



ft = 



PlVi 



where pi is the density of water at that pressure, and 
Vi the associated diver volume. Both p x and Vi are 
functions of pressure, and G will increase or decrease 
with increasing pressure depending upon whether the 
compressibility of water is greater than or less than 
that of the diver body. AG is not a perfectly linear 
function of AP. We have carefully measured the 
nonlinear term and found 

^,=8X 10- 8 g/cm Hg/cm Hg. 2 

Integrating twice gives us 

G= (4X 10- 8 ) P 2 +C 1 P+C 2 

where Oi and C 2 are the constants of integration in- 
troduced thereby. This term adequately represents 
the mass of water displaced by the diver as a function 
of pressure. 



2 Although, as previously stated, the pressure system has not been calibrated 
in absolute terms, within 1 percent, a load change of 1 kg on the separator plate 
results in a pressure change of 1 cm Hg. 



The experiment consists of loading the diver, first 
with two standards, Si and S2, and noting the equi- 
librium pressures, P S1 and P s2 associated with each. 
Following this, the unknown object, X, is placed on 
the diver and P x observed. The masses, M si and 
M S2 , and the volumes, y sl and V S2 , of the two stand- 
ards are assumed known. The equilibrium equation 
of the diver is 

M D ~G+M L - PW V L =0 

where M D is the diver mass, M L and V L the mass and 
volume, respectively, of the load and p w the density 
of water at the equilibrium pressure. As pointed 
out in the introduction M L —p w V L is the apparent 
mass of the load in water, B L . For each of the three 
loads placed on the diver during the experiment, we 
may write an equilibrium equation 

(A) Mz>-(4X lO-^FL-ftPsi-ft+Bs^O 

(B) M D -(AX10- 8 )P 2 S2 -C 1 P &2 -C 2 +B S2 =0 

(C) M D -(AX10^)P X -C 1 P X -C 2 +B X =0 

subtracting the second equilibrium equation from 
the first, and then the third from the first gives us 
two difference equations 

-(4X10- 8 )(Pl 1 -P| 2 )-ft(P sl -Ps2) + (5 sl -P S 2)=0 

-(4X10- 8 )(PI 1 -P^)-ft(P sl -P x ) + (P sl -P x )=0. 

The first difference equation may be immediately 
solved for a numerical value of d 



<7i= 



(p sl -p g2 )-4xio- 8 (Pi 1 -py 

(Psi~Ps2) 



and with this value of Oi inserted in the second 
difference equation, we may solve for the apparent 
mass of the unknown in water, B x . 

^= J B sl -4X10- 8 (PI 1 -Py-C 1 (P 81 -P^). 

This experiment demands a minimum of three 
observations — P S1 , P S2 , and P x . Our system has 
drift associated with electronic drift, drift in ambient 
conditions and (probably) other sources. In order 
to minimize the effect of this drift we use one of our 
loads as a drift monitor. This load is observed 
every second observation. For an experiment with 
only one unknown the observation format would 
contain five observations taken at uniform intervals, 
thus: P S1 — P S2 — Psi — Px — Pbi- The value of 
(Psi—Ps2), required in the calculations shown above 
is obtained by subtracting the observed value of P S2 
from the average of the first two P S i's. (Psi—Px) 
is obtained in a similar manner. The three ob- 
served values of P S i are averaged to give the effective 
value of P S i during the experiment. Effective 
values of P S2 and P x used in calculations of Oi and 
B xl are obtained by applying the difference terms 
to this effective value. 
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The values of water density used in the calculation 
of apparent mass in water of the standards (from 
assigned values of mass and volume) are based upon 
an assumed density of 0.999973 g/cm 3 at 4 °C and al 1 
atm, and corrected for departure from 1 at m of the effec- 
tive values of observed pressure. Errors in the assumed 
value of water density exert only a second-order 
effect on the calculated density of the unknown, 
because our density knowledge is primarily obtained 
from the density standards used in the experiments. 
The water is merely the comparison medium. 

The diver system is currently being used in several 
density studies not pertinent to this report. To 
illustrate typical diver operation, some data from 
these studies will be presented here. Table 1 shows 
the observed raw data taken during measurements 
on two imperfect silicon crystals, Xl and X2. The 
standards used in this experiment were two perfect 
silicon crystals, Si and S2, whose assigned density 
value was 2.3290040 g/cm 3 . The observed values of 
P E and the pressure differences were calculated from 
values of rising and falling pressures and velocities 
as previously explained. This table shows data 
taken during four independent experiments per- 
formed during January 1965, and the average value 
of a group of nonindependent similar experiments 
performed on the same crystals during the previous 
November. Between each experiment all crystals 
were removed from the hydrostatic system and dried 
and recleaned. 

Table 1 
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Table 2 shows the calculated values of effective 
P E and the assumed p w associated with these values. 



The apparent masses in water of the standards is based 
upon assigned true mass values of M S i = 2. 0000912 g 
and M S2 = 2. 0004884 g. The effective values of 
( i were obtained as previously described and they 
led to the values of apparent masses in water of the 
unknown cyrstals. 

Table 2 
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As pointed out in the introduction, we may 
obtain values of M x and V x from a simultaneous 
solution of the air weighing and the water weighing 
equations. Alternatively we may enter the water 
weighing equation with the value of B x obtained 
with the diver and a previously determined M x 
and solve for the density of the unknown, p x . 



Px~- 



Pw 



B x /M x -l 



where p w is the assumed water density at pressure P x . 
Very careful analysis of this point will show that the 
two methods are nearly identical. The values of 
p x calculated from the B x 's determined with the 
diver and M X1 = 2.0004774 g and M x2 =2. 0003478 g 
were: 



G January. 

7 January. 

8 January. 

9 January . 
Nov. Avg_ 

Average 



Crystal XI 



2. 3289950 
2. 3289958 

2. 3289944 
2. 3289950 
2. 3289956 

2. 3289951 



Crystal X2 



2. 3290012 
2. 3290018 
2. 3290017 
2. 3290015 
2. 3290019 

2. 3290016 



Pooled standard deviation = 4. 4X 10 _7 ~0.22 ppm. 
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5. Use of the Diver System — Density 
Changes 

There is widespread scientific interest in the 
measurement of slight density changes and density 
differences as well as in the measurement of absolute 
density itself. The Cartesian diver system is partic- 
ularly well-suited to determination of changes and 
differences because of its ability to provide high 
precision data on samples of small size. In problems 
such as determination of density differences be- 
tween nominally identical materials but having 
slightly different physical properties (such as be- 
tween hard drawn and annealed wire), or determina- 
tion of small density changes when a material is 
exposed to a density-modifying process (such as UV 
irradiation), homogeneity of the sample is usually 
of prime importance. Homogeneity can frequently 
be verified or inhomogeneity evaluated in samples of 
a few hundred milligrams, but it is frequently im- 
practical or impossible to do so in samples of several 
grams. Some density modifying processes can pro- 
vide homogeneous density changes only in small 
samples, as, for example, processes depending upon 
x-ray beams that have limited diameters or pene- 
trating powers. 



In measuring density changes in an object sub- 
jected to a density modifying process, the diver 
system is largely free of the requirements for high 
accuracy standards. If the two standards are cut 
from the same piece of parent material from which 
the object of interest is cut, and to about the same 
geometry, then we are justified in assuming that the 
handbook value of the density of the material applies 
to the standards absolutely. 

The diver was recently called upon to measure the 
change in density of a single crystal sample of rutile 
weighing less than 260 mg, when the stoichiometry 
was altered by thermal reduction. The standards 
used were cut from adjacent positions in the same 
boule from which the object of interest was cut. 
All three crystals were carefully weighed and two 
independent measurements on the apparent mass in 
water of the unknown, B x , were made based on the 
assumption that the density of the standards was 
4.249 g/cm 3 . 

The crystal of interest was partially reduced as 
called for in the experimental plan and was then re- 
weighed. Two additional observations were made 
upon B x . The data of interest is shown below: 



Weight of crystal of interest 

Two measured values of B X - 

Two calculated values of p X - 

Average p x 



Before modification 



0. 2586061 g 



After modification 



0. 2585972 g 



0. 1977429 
4. 248973 



0. 197426 
4. 248952 



0. 1977389 
4. 249169 



0. 1977388 
4. 249162 



4. 248962 



4. 249166 



Change of p x during modification +0. 000204 



The above data show the agreement in B x to be 
about 0.3 Mg, which is somewhat poorer than average 
diver operation. Nevertheless, although the change 
in density associated with the process under exami- 
nation amounted to only 5:10 5 , the system reported 
this change to agreement of about 7>2 percent. 

This experiment was performed upon a crystal 
of only 250 mg. With larger samples, the reproduci- 
bility of B x would remain about constant, hence the 
reproducibility of values of density would be much 
improved. 

6. General Comments 

Although our work was not primarily dhected to 
evaluation of diver performance as a function of 
diver characteristics, the following comments may 
be of interest: 

(A) We were never successful in working with 
Pyrex divers. Although the ones we tried out were 
always thoroughly annealed after fabrication, we 



never were able to achieve stability of much better 
than ±10 /ig. 

(B) The sensitivity of the sytem may be crudely 
denned as ( — )AB/AP e , that is the indicated change 
in apparent mass per unit change in equilibrium 
pressure. Defined thus, the sensitivity of the diver 
reported on was about ){i /xg/g change in P E . We 
experimented upon highly compressible forms in 
which sensitivity was 10 Mg/g? and also upon struc- 
tures with compressibilities so small that sensitivity 
approached 0.003 Mg/g- We found that the re- 
producibility decreased as the divers had extremely 
small or extremely large values of sensitivity. For 
the extremely sensitive structures we experienced 
difficulty in obtaining hydrostatic balance, and for 
the relatively insensitive structures, zero drift was 
most annoying. 

(C) The system is quite sensitive to vibration and 
the mechanical shock associated with closing doors, 
etc. Spaepen also reported this vibration sensitivity 
and described his efforts to minimize it. Our experi- 
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ment was conducted on a massive pier which had 
previously been used to support a high sensitivity 
balance used in the NBS mass calibration program. 
A vibration survey conducted prior to the experiment 
indicated that at 30 Hz, the vibrations were about 
6X10~ 4 g, decreasing somewhat at night. Our best 
work was done at night and on weekends. 

(D) Many experimenters working will) other 
buoyant devices have reported difficulty in achieving 
smooth "lift-off" of the buoyant structure from an 
alining cone at the bottom. We found that a blunt 
tip at the bottom of the quartz diver moved nicely 
out of a chromium-plated brass cone whose radius 
was somewhat greater than the quartz. We experi- 
mented with cones of both Pyrex and sapphire with 
poor results. 

(E) In a private communication Spaepen informed 
us that a newly fabricated diver is rarely stable — 
that a one or two months aging period is called for. 
This has also been our experience. We believe that 
annealing the quartz diver would probably reduce 
this aging period. 

(F) A major uncertainty in all hydrostatic weigh- 
ing is caused by the adherence of gas to the surfaces 
of the load. An air bubble large enough to displace 
1 jdg of water is easily detected visually and easily 
removed. An equal amount of gas dispersed among 
several hundred (or thousand?) bubbles would be 
impossible to detect (without a priori knowledge of 
the density of the material being hydros} at ieally 
weighed). This could lead to serious systematic 
errors in the values of density calculated from ob- 
served data. In order to randomize the surface 
trapped air, we always remove all loads (both 
standards and unknowns) between each hydrostatic 
experiment and thoroughly dry and reclean them. 
A stable volume of gas adhering to the surfaces of the 
diver body is less annoying, since ii> chief effect is to 
change the values of (\ and G which are determined 
in each experiment. 

(G) In changing diver loads the portion of the 
tools used which come into contact with the loads 
are kept continuously immersed in water, never 
allowing them to come into contact with air. In this 
way we minimize the likelihood of transferring air 
from the tools to the loads. 

(H) Generally, a perfectly cylindrical diver body 
will be insufficiently compressible for practical diver 



operation. Our method of achieving the desired 
compressibility is to measure AG/AP on a newly con- 
structed diver in its initial cylindrical form. We then 
partially flat ten the cylindrical walls by flame treat- 
ment in small increments, measuring AG/AP after 
each operation. This process is continued until the 
desired sensitivity is attained. Figure 1 shows this 
partially flattened section. 
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